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ABSTRACT 

Measurements of the temperature and density structure of the solar corona 
provide critical constraints on theories of coronal heating. Unfortunately, the 
complexity of the solar atmosphere, observational uncertainties, and the limi- 
tations of current atomic calculations, particularly those for Fe, all conspire to 
make this task very difficult. A critical assessment of plasma diagnostics in the 
corona is essential to making progress on the coronal heating problem. In this 
paper we present an analysis of temperature and density measurements above the 
limb in the quiet corona using new observations from the EUV Imaging Spec- 
trometer (EIS) on Hinode. By comparing the Si and Fe emission observed with 
EIS we are able to identify emission lines that yield consistent emission measure 
distributions. With these data we find that the distribution of temperatures in 
the quiet corona above the limb is strongly peaked near 1 MK, consistent with 
previous studies. We also find, however, that there is a tail in the emission mea- 
sure distribution that extends to higher temperatures. EIS density measurements 
from several density sensitive line ratios are found to be generally consistent with 
each other and with previous measurements in the quiet corona. Our analysis, 
however, also indicates that a significant fraction of the weaker emission lines ob- 
served in the EIS wavelength ranges cannot be understood with current atomic 
data. 

Subject headings: Sun: corona 



1. Introduction 

The origin of the high temperature plasma that permeates the solar corona has defied 
understanding for many decades. In principal, the thermal structure of the solar corona holds 
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many clues to the physical processes that convert magnetic energy into thermal energy. For 
example, it has been proposed th at the corona is he ated by frequent bursts of magnetic 



reconnection called nanoflares (e.g. |Parkerlll972l . 1 19831 ). In this model turbulent motions in 



the solar photosphere lead to the constant tangling and braiding of the magnetic fields that 
rise up into the corona. The dissipation of this topological complexity leads to the release 
of energy on very small spatial scales. This suggests that the corona should be composed of 
many fine loops that are in various stages of heating and cooling. Thus the distribution of 
temperatures and densities is a critical constraint on the frequency, duration, and magnitude 
of the heating events that give rise to the high temperature corona. 

Unfortunately, determining the distribution of temperatures and densities in the corona 
is a nontrivial problem. The solar corona is highly structured and highly dynamic, making 
it difficult to isolate individual structures. Obtaining accurate atomic calculations is also 
a problem. Because of its relatively high elemental abundance, emission from Fe has been 
the focus of many recent solar instruments. Interpreting observations from imaging instru- 
ments such as SXT /Yohkoh, EIT/SOHO, TRACE, EUVI/ STEREO, XRT /Hinode, and the 
upcoming AIA/ SDO depends critically on accurate atomic calculations for Fe. The accuracy 
of the available atomic calculations for this complex atom, however, is often unclear. 

The launch of the EUV Imaging Spectrometer (EIS) on Hinode has greatly expanded 
spectroscopic observations of the solar corona. EIS combines a broad temperature cover- 
age (Fevill-Fexvil, Fe xxn-Fe xxiv) with relatively high spatial (1") and spectral (22 mA) 
resolution, allowing the temperature and density structure of the corona to be examined in 
great detail. Of particular interest are the properties of coronal loops observed in solar active 
regions and flares. Current models of active region loops suggest that the observed densities 
and temperatures are signatures of non-equilibrium processes, and EIS, with its advanced 
diagnostic capabilities, provides much stricter observational constraints on physical models. 
Recent work on the properties of active region loops, for example, has suggested that ac- 
tive region loops near 1 MK h ave narrow distribu tions of temperature, high densities, and 



relatively small filling factors (IWarren et al.ll2008l ). 



Since many of the plasma diagnostics provided by EIS are based on Fe emission lines it is 
important to assess them critically and compare results with previous measurements. Some 
initial results have been alarming. For example, EIS spectroheliograms i n Fevill and Siv il 



are nearly identical, suggesting a very similar temperature of formation (jYoung et al.l 120071 ). 
The respective peaks in the ionization fractions, however, are separated by over 0.2 MK. 
This calls into question how accurately temperatures can be measured using the Fe emission 
observed in the EIS wavelength ranges. 



In this paper we present an analysis of EIS observations above the quiet solar limb. 
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Many previous observations of the quiet corona h ave suggested that the distribution of tem- 



peratures is very narrow, almost is o thermal fe.g.. iRavmond et al 
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2|). By comparing Fe and Si emission observed with EIS we are able to identify plasma 



diagnostics that are both self-consistent and in agreement with earlier results. Fortunately, 
most of the useful diagnostics that are identified are the strongest emission lines that can 
be observed with EIS. However, our analysis also shows that a significant fraction of the 
observed emission lines cannot be understood with current atomic models. 

In this work we also apply a Monte Carlo Markov Chain (MCMC) emission measure 
algorithm ( IKashyap fe Drakel Il998l . [2000) to the observed EIS spectra and find that the 
distribution of temperatures in the quiet corona is more complicated than previously thought. 
The emission measure is sharply peaked near 1 MK, as was found in earlier studies. We also 
find, however, a tail in the emission measure distribution that extends to higher temperatures. 
A high temperature component to the e mission measure distribut i on is a critical element of 
impu lsive coronal heating models (e.g., ICargill fe Klimchukl 120041 : iPatsourakos fe Klimchuk 



boosl ). 



Previous measurements of elec tron densities in the quiet corona above the limb have 



yielded values of logn e ~ 8.3 (e.g., iDoschek et al.lll997l ). EIS has several density sensitive 



line ratios that are useful in the quiet corona. We find that these ratios are all generally 
consistent with each other and with previous results, although there is considerable dispersion 
in the densities inferred from EIS. 

Finally, we also discuss the potential for measuring relative abundances in the corona 
using the Sx 264.233 A and Ovi 184.117 A lines observed with EIS. 



2. Instrumentation and Data Reduction 

The EIS instrument on Hinode produces stigmatic spectra in t wo wavelength r anges 



(171-212 A and 245-2 91 A) with a spectral resolution of 0.0223 A flCulhane et all 120071 : 



Korendyke et al.l 120061 ). There are 1" and 2" slits as well as 40" and 266" slots available. 
The slit-slot mechanism is 1024" long but a maximum of 512 pixels on the CCD can be read 
out at one time. Solar images can be made using one of the slots or by stepping one of the 
slits over a region of the Sun. Telemetry constraints generally limit the spatial and spectral 
coverage of an observation. For this work we focus on the results from a special observing 
sequence (HPW_FULLCCD_001) that returned the entire wavelength range of the CCD over a 
small region on the Sun (128" x 128"). At each position in the raster a 90s exposure is taken. 
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The data for this study were taken between August 19, 2007 23:47 UT and August 20, 
2007 03:01 at the west lim b. Context images from the EUV Imaging Telescope on SoHO 



( jDelaboudiniere et al.lll995l ) are shown in Figured] and indicate that the observed region was 
quiet during this time. Inspection of EIT 195 A movies during this period show no indication 
of significant coronal activity. The EIT 304 A image does indicate the presence of some cool 
material at low heights in the corona. 

These EIS data were processed using standard algorithms to remove the CCD pedestal, 
dark current, "cosmic ray" spikes, and warm pixels. The data numbers recorded in each pixel 
were also converted to physical units (erg cm -2 s _1 sr _1 A -1 ). For each intensity value an 
uncertainty is also calculated. This uncertainty includes counting statistics and read noise, 
but not the uncertainty associated with the absolute calibration. The implications of this 
will be discussed in later sections. 

There are several instrumental effects that impact our analysis. The first is the oscil- 
lation in the line centroids. This oscillation, which is about 2 spectral pixels in magnitude 
(~ 0.04A) over a period of about 90 minutes, is believed to be due to changing thermal 
conditions on the spacecraft during an orbit. To correct for this we assume that the average 
Doppler shift in the Fexn 195.119 A line averaged along the slit is zero at each slit posi- 
tion. Another important effect is the spatial offset between the two detectors. Because of 
a misalignment between the two CCDs there is a vertical offset of approximately 18 pixels 
between images taken in the different channels. In data taken before 2008 August 24 there 
is also an offset of 1-2 pixels in the solar-X direction. 

After these instrumental effects are accounted for we can determine the line intensity 
at each spatial position by calculating either moment of the line profile or a Gaussian fit. 
For making context rasters we use moments and EIS rasters in selected lines are shown in 
Figure [TJ These rasters also indicate the presence of some cool material at low heights in 
the corona. 

Since our objective is to analyze observations from the quiet corona with very high signal 
to noise we have computed a spectrum averaged over a 35" x 104" region that lies above 
this cool material. The center of this region is about 70" above the limb or about l.O7i? . 
In constructing this average spectrum we have ignored any pixel that has been marked as a 
warm pixel or as a cosmic ray impact. The intensity in each spectral pixel is computed from 
the average of approximately 3400 spatial pixels. 



We have used the identifications of iBrown et al.l (120081 ) and I Young et al.l (120071 ) to de 



termine the emission lines of interest. For each of these lines we have fit the line profiles 
with Gaussians and extracted the relevant line intensities. These line intensities are given 
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in Tables [T] and El Ad ditiona ll y, we extracted line intensities for three Fe IX lines that were 
recently identified by lYoungl (120081 ). This list is not complete. The objective here is to 
consider the strongest emission lines that could be used in spatially resolved observations, 
and not to consider every emission line in the spectrum. 

All of the lines except Fexvi 262.984 A are well represented by Gaussians. There is 
very little signal in the Fexvi profile and the observed emission is essentially noise. This 
measurement does provide a very useful upper bound on the amount of high temperature 
emission in this region. 

Tables [1] and [2] give the most significant atomic transitions for each emission line. The 
CHIANTI level numbers are also given in these tables. The level numbers are simply an 
ordered list of the transitions. They aid in identifying which emission lines involve transitions 
to the ground state as well as which lines originate in the same upper level. The Fex 190.038 
and 184.536 A lines, for example, originate in the same upper level and form a branching 
ratio. The branching ratio only depends on the relative decay rates and should be more 
accurate than other ratios. In these tables the emissivity at the peak of the ionization 
fraction and a density of logn e = 8.35 is also given. 

In addition to the Fe lines EIS also observes several weaker Sivil, IX, and x emission 
lines. Si emission lines from these ionizati on stages have been used in previous emission 



measure analysis of the quiet c orona (e.g., iFeldman et al J 1 19991 ; IWarrenlll999l : lLandi et al. 



20021 : IWarren fe Warshalll 12002 ) and provide a useful comparison for the analysis of the Fe 
lines. 

Elemental abundances play an important role in determining the magnitude of radiative 
losses in the corona. There are several emission lines from high first ionization potential 
elements observed within the EIS wavelength ranges (IFeldman et al.l 120081 ). At coronal 
temperatures the Sx 264.233 A and Sxill 256.686 A lines can be used for studying the 
composition. Only the S x line appears in these data. The O VI lines at 183.937 and 184.117 A 
provide another measurement from a high FIP element. Ovi is Li-like and the ionization 
fraction for this ion is significant at coronal temperatures. 



3. Temperature and Density Measurements 

The observed intensities are related to the plasma emissivities, e\{n ) T) ) and the differ- 
ential emission measure, £(T), by the expression 



47T 



e A (n e ,T)£(T) dT. 



(1) 
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In this context the plasma emissivities are the radiated power (erg s _1 ) divided by the square 
of the electron number density, n e . For many atomic transitions this quantity depends 
mainly on the electron temperature, T. For many of the Fe lines considered here, however, 
the emissivity is strongly dependent on the density, even for transitions to the ground state, 
and to facilitate intensity calculations we have computed grids of emissivities over a wide 
range of densities and temperatures using the CHIANTI 5.2.1 atomic physics database (e.g., 
Landi et al.ll2006l ). The 5.2 .1 version of the data base corrects an error in the atomic data for 
Fexm. The abundanc es of iFeldman et al.l (119921 ) and the low density ionization fractions of 
Mazzotta et al.l (119981 ) are assumed. In this expression the emission measure is the line of 
sight emission measure, £(X) = n 2 e ds/dT, and has units of cm -5 K _1 . In the plots we will 
generally display the differential emission measure times the temperature. 

For this work we will consider three different methods for reconstructing the differential 
emission measure from the EIS intensity measurements. The first two methods rely on a 
parametrization of the emission measure. The simplest approximation is that of a single 
temperature plasma where the differential emission measure is a delta function 



£(T) = EM 5(T-T ). 



(2) 



To account for the possibility that there is some dispersion in the temperature distribution 
we also consider a Gaussian representation of the differential emission measure 



EMn 



exp 



(T - T ) 
2o\ 



21 



(3) 



The Gaussian DEM is parametrized so that for very narrow temperature distributions we 
recover the parameters for the isothermal case, i.e., for T§jo T 3> 1 we have J £(X) dT ~ EM . 
To determine the best-fit parameters for either of these emission measure models we use a 
Levenberg-Marquardt technique implemented in the MPFIT package. We have implemented 
this algorithm so that the density in Equation [I] can either be a free parameter or have a 
fixed value. 

Fina lly, we also apply a Monte C arlo Markov Chain (MCMC) emission measure al- 
gorithm (IKashyap fe Drakd Il998l . I200G ) distributed with the PINTof ALE spectral analysis 
package to these data. This algorithm has the advantage of not assuming a shape for the 
differential emission measure. The MCMC algorithm also provides for estimates of the error 
in the DEM. In its current implementation the MCMC algorithm does not allow for the 
density to be a free parameter. 

Perhaps the most important test of the EIS spectra is the application of the isothermal 
DEM model to the observed Si emission. The atomic data for Si appears to be very well un- 
derstood and previous emission measure analysis in the quiet corona has yielded consistent 
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results (e.g.. iFeldman et al.lll999l : IWarren fe Warshallll2002l ). The isothermal DEM calcula- 
tion for the Si lines is illustrated in Figure El where the best-fit parameters are shown. Here 
we also display the emission measure loci curves defined by 



EM(T) = 



(4) 



Previous analysis in the quiet corona above the limb has shown that these curves tend to 
intersect at a point, suggesting isothermal plasma. As noted before, we solve for the best fit 
parameters thro ugh y 2 minimization rather than estimating them by eye as has been done 
previously (e.g.. IFeldman et al.lll998l ). 



The best fit temperature is similar to the values o f 



or Tp = 6.05 determined bv lLandi et al. 



(120021 ) and log T = 6.01 determined bv [Allen et al.l (12000r). Other studies yield somewhat 
higher temperatures of log T n ^ 6.15 (e.g., Raymond et al.l 1997 : Feldman et al. 1999 : Warren 



1999l : IWarren fc Warshallll2002l ) . Additionally, as indicated by the values listed in Tabled the 
computed intensities of the Si VII and Si x lines are consistent with each other to within about 
±15%. This calculation depends on the density, which we have assumed to be logn e = 8.35 
in this calculation. If we allow the density to be a free parameter in the minimization then 
we find logn e = 8.22. These values are within the range of possible densities. We will discuss 
the densities derived from the various line ratios and DEM inversion methods in detail at 
the end of this section. 

The generally good agreement among the observed Si lines is not matched by the ob- 
served Fe spectrum. Emission measure loci plots for all of the Fe lines, which are shown 
in Figure [3l do not reveal any discernible pattern and do not suggest either isothermal or 
multi-thermal plasma. 

If the atomic data, the assumed density, and the observed intensities were mutually 
consistent then all of the emission measure loci curves for a given ion would lie very close 
together, as they do for Sivil and Six. Closer inspection of the the emission measure loci 
curves for Fe shown in Figure [3] suggests that for each ion there are several emission lines 
that are mutually consistent and others that are discrepant by varying amounts. Here we 
identify these lines and discuss previously identified blends. 



Fevill The emission measure loci for 185.213 an d 186.601 A li n es are consistent, but the 
curve for 194.663 A is about a factor of 2 higher. iBrown et al.l (120081 ) indicate a possible 
blend of Fevill 194.663 A with O V 195.593 A, which would have negligible intensity at this 
height above the limb. In the core of an active region Fevill 185.213 A is blended with 
Nixvi 185.251 A and Fevill 186.601 A is blended with Caxiv 186.610 A. As stated earlier, 
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the ionization fraction for FeVlll may need revision, so these lines should be used with 
caution. 



Fe IX The newly identified 188.497, 189.941, and 197.862 A lines are all in good agreement. 
The emission measure loci for Feix 171.073 A is off by about a factor of 2. The effective area 
for EIS at this wavelength is very low, making the observed intensity highly uncertain. 



Fex The 174.532, 177.239, and 184.536 A lines are all in agreement. The Fex 190.038, 
207.449, and 257.262 lines are not. The 184.536 and 190.038 A lines originate in the same 
upper level and form a branching ratio. The theoretical ratio, however, is 3.56 while the 
observed ratio is 2.70. The Fex 257.262 line forms a density sensitive ratio with 184.536 
and 190.038 A, but the problems with the emission measure loci suggests that the densi- 
ties derived from this ratio are not consistent with the densities derived from the Si ratio. 



Brown et all (120081 ) indicate a possible blend of Fex 184.536 A with Arxi 184.524 A, which 



would be a problem for active region observations. 



Fexi The 180.401, 182.167, 188.216, and 192.813 A lines are all in agreement. The emission 
measure loci for 188.299, 257.547, and 257.772 A dif fer by factors of 2 to 4. The 182.167 and 
188.216 A lines form a density sensitive line ratio. iBrown et al.l ( 120081 ) indicate that Fexi 
180.401 A is blended with Fex 180.407 A. Using the isothermal emission measure determined 
from Si to estimate the intensities of the two lines suggests that the Fe x contribution to the 
observed intensity is about 10%. There are O V lines that can make important cont ributions 
to the observed Fexi 192.813 A line profile during transient events ( Iko et aL 20091 ) . 



Fexn The Fe xii 186.880, 192.394, 193.509, 195.119, and 196.640 A lines are all in relatively 
good agreement. The 203.720 and 256.925 A lines differ by factors of 2 and 30, respectively. 
The 186.880 and 196.640 A lines form density sensitive ratios when paired with any of the 
192.394, 193.509, and 195.119 A lines. T here is a S xi line at 186.84 A, but it is ge nerally weak 
compa red with the Fexn 186.880 line ( jYoung et al.l 120081 ). In active regions, I Young et al. 



( 120081 ) note that there is an Fexn 195.18 A line that becomes important at high densities. 
This component is generally small but can impact Doppler shift and line width measurements 
with the Fexn 195.119 A line. 



Fexm The 196.525, 197.434, 200.021, 202.044, and 203.826 A lines are generally consistent. 
The emission measure loci for 201.121, 204.937, 246.208, and 251.953 are all offset by varying 
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amounts. The 196.525 and 203.826 A lines form density sensitive ratios with 202.044 A. There 
is a possi ble blend of Fe xm 1 96.525 A with Fevil l 196.65 A, but the Fevill line is believed to 
be weak dYoung et alJl20Q8| ). brown et all d2008h indicate that Fexm 201.121 A is blended 
with Fexn 201.121 A. 



Fexiv The 211.316, 270.519, and 274.203 A emission measure loci curves are all relatively 
close to each other while the curve for 264.787 A is offset. The 264.787 and 274.203 A lines 
form a density sensitive pair, but this ratio should be in the l ow density limit in the quiet 
Sun. There is a blend of Fexiv 274.203 A with Sivil 274.175 A. [Young et all hoQlh indicate 
that the Si VII 274.175 A intensity is less than 0.25 times the intensity of Si VII 275.352 A. The 
ratio of the Si lines is sensitive to density, however, and for the the low densities measured 
above the limb the ratio is calculated to be about 10%, suggesting a negligible contribution 
of Si VII to the Fexiv 274.203 A intensity measured here. 



Fe XV and Fe xvi There is only a single Fexv emission line in the EIS wavelength ranges. 
There are several Fexvi lines present, but as indicated earlier, no emission is observed in any 
of these lines. The intensity given in Tabled] for Fexvi 262.984 A serves only as an upper 
bound and limits the peak temperature in the emission measure. There is no independent 
verification of the consistency of these lines. 

We have used the lines identified here by their internal consistency to compute the 
emission measure using the three different models. Both the emission measures and the 
emission measure loci are shown in Figure BJ The isothermal solution for these selected Fe 
lines is very similar the solution found for the Si lines. The best fit temperature is logT = 
6.05, compared with 6.07 for the Si lines. It is clear, however, that the isothermal model 
cannot reproduce both the low temperature emission (Fe ix-xiii) and the emission observed 
in the higher temperature lines (Fexiv-xvi). The intensities for the high temperatures lines 
computed from the isothermal emission measure are systematically too small by a factor of 
about 4. 

This systematic discrepancy suggests that the plasma in the quiet corona is not isother- 
mal. To investigate this we have applied the Gaussian DEM algorithm to the Fe lines. If 
we include only the Feix-xill lines we obtain a very narrow emission measure distribution 
(gt — 4.8). Including the Fexiv-xvi lines leads to a somewhat broader emission measure 
distribution (a T = 5.2), but only a marginal improvement between the observed and calcu- 
lated intensities. For this case the observed intensities are systematically larger than what 
is calculated from the Gaussian emission measure by about a factor of 2. 
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The failure of the Gaussian emission measure model to reproduce all of the observed 
line emission is due in large part to the symmetry that this model imposes on the DEM. 
Increases in the emission measure at higher temperature must be accompanied by increases 
at lower temperatures, and this would lead to discrepancies in the Feix and Fex intensities. 
The MCMC algorithm does not assume any functional form for the emission measure and 
provides for more flexibility. The differential emission measure calculated using the MCMC 
algorithm is displayed in Figure [4] and calculated intensities are given in Table 0J In this 
case the agreement between the observed and calculated intensities is improved for the high 
temperature lines while maintaining the good agreement at the lower temperatures. This is 
achieved by the introduction of a high temperature tail in the differential emission measure 
in addition to the strong peak near 1 MK. The strong peak in the emission measure is 
consistent with previous work. Previous studies have usually employed relatively restrictive 
parametrizations for the emission measure and would not have been able to identify this 
high temperature component. 

The best fit parameters for the Si and Fe isothermal emission measure calculations 
show a significant discrepancy in the magnitude of the emission measure, with the Si emission 
measure being about 60% larger (logE'Mo = 26.49 for Fe and 26.72 for Si). We find a similar 
difference if we use the MCMC DEM to infer the intensities of the Si lines. Curiously, the Si 
lines all fall on the long wavelength detector while the Fe lines all fall on the short wavelength 
detector. None of the long wavelength Fex, XI, XII, or xill lines are consistent with the short 
wavelength lines from the same ion. Similarly, the high temperature lines, whose intensities 
the isothermal models systematically under-predict, generally lie on the long wavelength 
detector. This suggests that calibration differences between the two detectors might explain 
at least part of these discrepancies. While this explanation may sound plausible it does 
not withstand closer scrutiny. Perhaps most importantly, the magnitude of the discrepancy 
between the calculated and observed Fe xiv-xvi intensities is much larger than the difference 
between the Si and Fe isothermal emission measures. If the observed values were reduced by 
a factor of 1.6 the calculated intensities would still be off by a factor of 3 or more. One Fexiv 
line does appear on the short wavelength detector (211.316 A) and its behavior is similar to 
the long wavelength Fexiv lines. Finally, the discrepancies in the emission measure loci 
curves discussed previously do not show any systematic variation with wavelength. With 
these data we find no evidence that there is a problem with the relative calibration of the 
detectors. 

When computing the line intensity the density and temperature both play an im- 
portant role in determining the emissivity. This makes it difficult to determine the tem- 
perature and density independently. One way of addressing this issue is to allow both 
parameters to vary while using Equation [1] to find the optimal parameters. We have 
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done this in applying the isothermal model to the Si and Feix-xill lines and we gener- 
ally find somewhat different densities. The EIS wavelength ranges contain a large num- 
ber of density sensitive line ratios (see I Young et al.l 120071 for a summary) and we can 
use these data to compare the results. In Figure [5] we show the densities derived di- 
rectly from Six 258.375/261.058 A, Fex 184.536/257.262 A, Fexi 182.167/188.216 A, Fexn 
186.880/195.119 A, Fexn 196.640/195.119 A, Fexm 196.525/202.044 A, Fexm 203.826/202.044 A, 
and Fexiv 264.787/274.203 A. The emissivities are computed assuming logT = 6.05, as de- 
rived from the isothermal emission measure analysis. 

These ratios show some significant differences. The Fex and Fexiv ratios give useless 
results, with the Fex being below the low density limit and Fexiv yielding a density an order 
of magnitude larger than the others. The other ratios yield values ranging from logn e = 8.14 
to 8.50. These values are generally consistent with t he densities determine d in earlier analysis 
of quiet Sun limb observations, e.g., logn e <~ 8.3 by lDoschek et al.lll997l and logn e ~ 8.5 by 
Landi et al J 12002 . Based on previous analysis, the Six ratio is probably the most reliable 
ratio of this group. We note that the densities from the Fe lines are scattered around the 
density derived from Si and if we average the Fe densities together we obtain 8.34 ± 0.16. 
This is a variation of about a factor of 2. To simplify the discussion, all of the plots and 
tables assume a density of logn e = 8.35. 

The large region that we have used for compu ting an average sp ectrum leads to very 
small statistical errors for the measured intensities. lLang et al.l (120061 ) calculate the uncer- 
tainty in the absolute radiometric calibration for EIS to be 22%. This additional uncertainty 
can be easily added in quadrature to the values given in Tables [1] and [3 We have rerun 
all of the dem calculations including the calibration uncertainties and obtained very simi- 
lar results. The biggest differences are for the MCMC algorithm, which yields a somewhat 
smoother emission measure at high temperatures. Variations in the detailed structure of 
the emission measure a re not surprising given the ill-posed nature of the emission measure 
inversion problem (e.g. JCraig fe Brownlll976l ). 



Abundances 



To accurately compute the magnitude of the radiative losses in the corona we must have 
some measure of the coronal composition. Unde rstanding radiative c ooling is important for 



mode ling the evolution of coronal loops (e.g., iWarren et al.l 120031 : IWinebarger fe Warren 
20041 ) . Spatially averaged measurements suggest that the low first ionization potential ele- 
ments, such as Fe, Si, and Mg, are all enriched in th e corona relative to high FIP elements, 
such as C, N, and O by about a factor of 4 (e.g., iFeldman et al.lll998l ). Temporally re- 
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solved meas urements in active region s indicate, however, that this enhancement may change 
with time (jWiding fe Feldmanll200ll ). with newly emerged active regions having more pho- 
tospheric abundances (i.e., no enrichment of the low FIP elements). This suggests the need 
for more systematic abundance measurements. 

Several abundan ce diagnostics in the EIS wavelength ranges have been discussed by 
Feldman et al.l (120081 ). For coronal plasma the most useful EIS emission lines from high 
FIP elements are Sx 264.233 A and Sxm 256.686 A. In these quiet Sun data only the Sx 
264.233 A line appears. The FIP for S is 10.3 eV and it lies at the bounda ry between the 
high and low FIP elements. Some measurements (e.g., IFeldman et al.lll998r) sh ow a partial 



Feldman et al 



19921 ) indicate es- 



enrichment of S in the corona. Other measurements (e.g., 
sentially photospheric abundances. The general trend is for S to behave somewhat differently 
than the low FIP elements, and S x can be used as a proxy for measuring relative abundance 
variations in EIS observations of the corona. For these observations above the quiet corona 
it is also possible to use O VI as a high FIP line (e.g., iFeldman et al.lll998l ). O VI is Li-like 
and the ionization fraction extends to high temperatures. 

In calculating the emissivities we have assumed a coron al composition with A Fe = 8.10, 
A .Qj = 8.10, An = 8.89, and A$ = 7.27 (IFeldman et al.lll992l ). The photospheric composition 
of lGrevesse fe Sauvall dl998h has A Fe = 7.50, A Si = 7.55, A Q = 8.83, and A s = 7.33. If we 



use the MCMC emission measure computed from the Fe lines to compute the the intensities of 
O VI 183.937, 184.117 and Sx 264.233 A we obtain values of 0.8, 1.6 and 15.0 erg cm" 2 s" 1 sr" 1 
respectively. The observed intensities in this re gion are 2.8, 4.8 and 34.8 erg cm -2 s _1 sr _1 , 
suggesting that the relative abundances of the (IFeldman et al.lll992l ) are not consistent with 
these data. The abundance of S and O in the corona would need to be increased by about 
a factor of 2.9 or the abundance of Fe decreased by about a factor of 2.9 to resolve this 
discrepancy. Similarly, the difference in the emission measures derived from the Si and Fe 
lines could be resolved by adjusting their relative abundances. 



5. Quiet Sun Disk 

The observations taken above the limb provide an opportunity to study a high signal to 
noise spectrum derived from an average over a largely homogeneous region. This analysis has 
revealed many problems with interpreting the observed intensities. These problems could 
be due to blends with other lines, errors in the calibration, or blends with other lines. Since 



the m orphology of the quiet solar atmosphere changes with temperature (e.g.. IFeldman et al 
Il999l ) examining spatially resolved images in these emission lines offers a means of identifying 
previously unidentified blends in the lines. 
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Most EIS observations preserve the information from only a few narrow spectral win- 
dows. Observations that contain the full wavelength range are typically small in size. We 
have found several observations that contain the full wavelength range and cover a relatively 
large region on the Sun (up to 128" x 512"). For one of these observations (taken November 
15, 2007 11:14 UT) we have processed the data and constructed rasters in as many of the 
emission lines given in Tables [1] and [2] as possible. These rasters are computed by performing 
Gaussian fits to the spectral lines at each spatial position. Unfortunately, some emission 
lines from the limb spectra are too weak to fit in the spatially resolved disk spectra. 

The rasters for many of the emission lines of interest are shown in Figures [6] and 
These rasters suggest a progression from small scale structures that correspond to the mag- 
netic network at low temperatures to much longer loops, with no network pattern at high 
temperatures. Surprisingly, there is little evidence for blending in many of the problem emis- 
sion lines. For example, the Fex 184.536 and 190.038 A rasters are nearly identical. The 
contrast in the Fex 257.262 A raster is somewhat different than in the other rasters, but 
this line is sensitive to density and this is to be expected. Similarly, the Fexi 188.216 and 
188.299 A rasters are nearly identical. There is some evidence for small bright features in 
the Fexi 192.813 A raster, suggestive of Ov emission. For Fexn the discrepant 203.720 A 
rasters appears to be consistent with the rasters from the other lines. The Fexn 256.925 A 
line does appear to be blended with a cooler line. The contrast between the center of this 
raster and the structures in the north and south is smaller than for the other Fexil lines. 
For similar reasons, the Fexiv 264.787 A also appears to be blended with a cooler line. 



Discussion 



We have presented a detailed analysis of spectroscopic observations taken in the quiet 
corona above the limb with the EIS instrument on Hinode. The strongly peaked differential 
emission measure calculated from the Si and Fe lines are generally consistent with each 
other and with previous measurements. This suggests that despite the complexity of the 
Fe atom, accurate differential emission measure calculations are possible with the strongest 
emission lines observed with EIS. Our ana lysis validates some of the initial EIS e mission 
measure calculations for coronal loops, (e.g.. [Warren et al.ll2008l : iTripathi et al.ll2009l ). which 
emphasize the lines given in Table @J 

Using the MCMC DEM algorithm we find evidence for a tail in the emission measure 
that extends to high temperatures. Such features are extremel y important to theories of 



coronal heating; based on the impul sive release of energy (e.g. ICargill fe Klimchuk 



2004 



Patsourakos fe Klimchukl 120061 . 120081 ) . It will be interesting to see if similar features of the 
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DEM are observed in active region loops. 

The densities measured from 6 of the 8 density sensitive line ratios that we considered 
are generally consistent with each other and with previous measurements. There is, however, 
considerable scatter in the measured densities. Measurements in active regions suggest that 
these differences are systematic, with the Fexn 186.880 /195.119 A density being systemati- 



cally higher than the Fexm 203.826/202.044 A density flYoung et all [20071 ). 



This work has highlighted problems with the analysis of many emission lines in the 
EIS spectral wavelengths. These differences result from some combination of errors in the 
atomic data, blends with unidentified lines, and uncertainties with the EIS calibration. Some 
of these difficulties are clearly due to the errors in the atomic data or assumptions that have 
been made in using them. The observations suggest that the ionization fractions for Fevill, 
for example, may need to be shifted to high temperatures. The quiet Sun disk rasters for 
Fex 190.038 A and 257.262 A, Fexi 188.299 A, Fexn 203.720 A, and Fexm 251.953 A do not 
appear to be blended with other lines. These rasters do suggest blends for Fexn 256.925 A 
and Fexiv 264.787 A. The possibility of errors in the EIS relative calibration is unclear. The 
Fexiv 211.316 A line is consistent with Fexiv 270.519 and 274.203 A, suggesting that there 
are no significant discrepancies. It is unsettling, however, that none of the Fex, XI, XII, and 
xill lines near 250 A agree with the lines at shorter wavelengths. 

The use of EIS observations for abundance measurements is still unsettled. The ratios of 
the S and O lines to Si and Fe do not yield results that are consistent with existing sets of solar 
abundances. It is unclear how these issues can be resolved at present. Despite this difficulty 
it is still possible to use EIS observations to investigate spatial and temporal variations in 
the solar composition. This is illustrated in Figure [6l where we can see spatially resolved 
images in Sx and Six. When combined with emission measure analysis, observations such 
as these can be used to understand the how much the composition varies from structure to 
structure in the solar corona. Work on the analysi s of quiet Sun disk spectra are currently 



in progress and will be reported in a future paper Brooks et all (I2009h . 



Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
domestic partner and NASA and STFC (UK) as international partners. It is operated by 
these agencies in co-operation with ESA and NSC (Norway). 
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Fig. 1. — (top panels) EIT images taken during the EIS full CCD observations. The boxes 
indicate the EIS field of view and the region used for computing the average spectra, (bottom 
panels) EIS spectroheliograms in selected emission lines. The dark vertical band in each 
image is due to atmospheric absorption during a brief orbital eclipse. 
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Fig. 2. — Emission measure analysis of the Si emission lines observed with EIS. In these 
calculations the density is held constant at logn e = 8.35. 
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Fig. 3. — Emission measure loci plots for Fevm xvi emission lines observed with EIS. 
These plots illustrate the problems with a number of the emission lines observed within this 
wavelength range. The red dot is the isothermal emission measure derived from the Si lines. 
In these calculations the density is held constant at logn e = 8.35. 
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Fig. 4. — Differential emission measures derived from Feix-xvi and different emission mea- 
sure models (isothermal, Gaussian, and Monte Carlo Markov chain). For the Gaussian and 
MCMC models the temperature times the differential emission measure is displayed. The 
MCMC DEM has the lowest x 2 and best reproduces the observed intensities. In these 
calculations the density is held constant at logn e = 8.35. 
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Fig. 5. — Electron densities derived from three density sensitive line ratios. The emissivities 
were evaluated at logT = 6.05. 



- 22 - 




Fig. 6. — EIS rasters from observations of the quiet solar disk taken November 15, 2007 
from 11:13:59 to 14:12:17 UT. The 1" slit has been stepped over a region 128" x 512" in size, 
and an 80s exposure has been taken at each slit position. Rasters for O V - Six are shown. 
These spatially resolved disk measurements allow for the morphology of different emission 
lines from the same ion to be compared. The Fexi 192.813 A shows a small contribution 
from Ov. 
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Fig. 7. — The same region as is shown in Figure El except that Fexn - Fexv are shown. The 
Fexii 256.925 A and Fexiv 264.787 A lines appear to have excess emission in the center of 
the rasters, suggesting blends with cooler emission lines. 



-24- 



Table 1. Fe Intensities Measured in the Quiet Corona above the Limb with EIS ; 
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25.02 


± 


0.16 










3s 2 3p 2 3 P 2 




3s 2 3p 3d 3 D 2 


3 


25 


1.35e-24 




± 






Fe xiii 


204.937 


6.19 


3s 2 3p 2 3 P 2 




3s 2 3p 3d 3 Di 


3 


23 


1.13e-24 


8.19 


± 


0.17 




Fe xiii 


246.208 


6.19 


3s 2 3p 2 3 Pi 




3s 3p 3 3 Si 


2 


14 


2.35e-24 


9.36 


± 


0.18 




Fe xiii 


251.953 


6.19 


3s 2 3p 2 3 P 2 




3s 3p 3 3 Si 


3 


14 


4.52e-24 


19.33 


± 


0.24 




Fexiv 


211.316 


6.27 


3s 2 3p 2 P 1/2 




3s 2 3d 2 D 3/2 


1 


11 


1.00e-23 


39.47 


± 


0.51 




Fexiv 


264.787 


6.27 


3s 2 3p 2 P 3/2 




3s 3p 2 2 P 3/2 


2 


10 


3.72e-24 


20.72 


± 


0.08 




Fexiv 


270.519 


6.27 


3s 2 3p 2 P 3/2 




3s 3p 2 2 P 1/2 


2 


9 


2.56e-24 


6.96 


± 


0.06 




Fexiv 


274.203 


6.27 


3s 2 3p 2 P 1/2 




3s 3p 2 2 S 1/2 


1 


8 


5.54e-24 


18.31 


± 


0.08 




Fexv 


284.160 


6.33 


3s 2 % 




3s 3p x Pi 


1 


5 


2.88e-23 


21.20 


± 


0.14 


? 
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Table 1 — Continued 



Ion 


A 


Tmax 


Transition LI 


L2 


eA 


lobs 


&Iobs 


Fexvi 


262.984 


6.43 


3p 2 P 3 /2 3d 2 D 5/2 3 


5 


9.18e-25 


0.42 ± 


0.04 ? 



a In this and all subsequent tables wavelengths are in A and the units for the intensities and corre- 
sponding uncertainties are erg cm -2 s _1 sr _1 . Lines that can be used for emission measure analysis are 
indicated by the symbol. There is no independent check on the Fexv 284.160 A and Fexvi 262.984 A 
lines so they are marked with a "?" . 



- 26 - 



Table 2. Si, S, and O Intensities Measured in the Quiet Corona above the Limb with EIS 



Ion 


A 


Tmax 


Transition 




LI 


L2 




lobs 




°Iobs 




Ovi 


183.937 


5.47 


Is 2 2p 2 P 1/2 


Is 2 3s 


2 s 1/2 


2 


4 


4.66e-25 


2.81 


± 


0.05 


? 


Ovi 


184.117 


5.47 


Is 2 2p 2 P 3/2 


ls 2 3s 2 S 1/2 


3 


4 


9.36e-25 


4.82 


± 


0.10 


? 


Si vii 


275.352 


5.77 


2s 2 2p 4 3 P 2 


2s 2p 5 


3 P 2 


1 


6 


4.55e-24 


11.54 


± 


0.07 




Si vii 


275.665 


5.77 


2s 2 2p 4 3 Pi 


2s 2p 5 


3 Pi 


2 


7 


7.04e-25 


1.81 


± 


0.05 




Si vii 


278.445 


5.77 


2s 2 2p 4 3 Pi 


2s 2p 5 


3 P 2 


2 


6 


1.45e-24 


2.95 


± 


0.19 




Si ix 


258.073 


6.05 


2s 2 2p 2 X D 2 


2s 2p 3 


X D 2 


4 


13 


1.31e-25 


5.08 


± 


0.08 




Six 


253.791 


6.13 


2s 2 2p 2 P 1/2 


2s 2p 2 


2 P 3 /2 


1 


10 


4.21e-25 


12.84 


± 


0.13 




Six 


258.375 


6.13 


2s 2 2p 2 P 3/2 


2s 2p 2 


2 P 3 /2 


2 


10 


2.19e-24 


71.30 


± 


0.18 




Six 


261.058 


6.13 


2s 2 2p 2 P 3/2 


2s 2p 2 


2 Pl/2 


2 


9 


1.30e-24 


43.34 


± 


0.13 




Six 


271.990 


6.13 


2s 2 2p 2 P 1/2 


2s 2p 2 


2 s 1/2 


1 


8 


9.76e-25 


37.62 


± 


0.10 




Six 


277.255 


6.13 


2s 2 2p 2 P 3/2 


2s 2p 2 


2 s 1/2 


2 


8 


7.99e-25 


28.83 


± 


0.11 


<? 


Sx 


264.233 


6.15 


2s 2 2p 3 4 S 3/2 


2s 2p 4 


4 P 5 /2 


1 


6 


8.21e-25 


34.77 


± 


0.11 


? 
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Table 3. Isothermal Model Applied to the Si Lines 



Ion 


Wavelength 


Icalc 


lobs 


Icalc/ lobs 


Si VII 


275.352 


11.41 


11.54 


0.99 


Si VII 


275.665 


1.83 


1.81 


1.01 


Si ix 


258.073 


5.71 


5.08 


1.12 


Six 


253.791 


14.32 


12.84 


1.12 


Six 


258.375 


74.43 


71.30 


1.04 


Six 


261.058 


45.39 


43.34 


1.05 


Six 


271.990 


34.09 


37.62 


0.91 


Six 


277.255 


27.89 


28.83 


0.97 
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Table 4. MCMC Model Applied to the Feix-xvi Lines 



Ion Wavelength I ca i c I obs I ca ic/Iobs 



Feix 


188. 


.497 


34. 


.13 


31. 


.28 


1. 


.09 


Feix 


189. 


.941 


18. 


.02 


15. 


.36 


1. 


.17 


Feix 


197. 


.862 


24. 


.50 


21. 


.02 


1. 


.17 


Fex 


174. 


.532 


560. 


.17 


572, 


.57 


0. 


.98 


Fex 


177. 


.239 


304. 


.81 


308. 


.28 


0. 


.99 


Fex 


184. 


.536 


119. 


.94 


142, 


.17 


0. 


.84 


Fexi 


180. 


.401 


415. 


.93 


432 


.24 


0. 


.96 


Fexi 


182. 


.167 


53. 


.84 


58. 


.50 


0. 


.92 


Fexi 


188. 


.216 


196. 


.01 


224 


.90 


0. 


.87 


Fexi 


192. 


.813 


50. 


.16 


57, 


.73 


0. 


.87 


Fexn 


186. 


.880 


33. 


.11 


35 


.17 


0. 


.94 


Fexn 


192. 


.394 


88. 


.03 


79, 


.48 


1. 


.11 


Fexn 


193. 


.509 


185. 


.60 


177, 


.53 


1. 


.05 


Fexn 


195. 


.119 


277, 


.87 


274, 


.67 


1. 


.01 


Fexn 


196. 


.640 


9. 


.04 


11, 


.03 


0. 


.82 


Fe xm 


196. 


.525 


1. 


.71 


2. 


.71 


0. 


.63 


Fe xm 


197. 


.434 


5. 


.24 


7. 


.10 


0. 


.74 


Fe xm 


200. 


.021 


6. 


.37 


9. 


.43 


0. 


.68 


Fe xm 


202. 


.044 


126. 


.16 


157, 


.66 


0. 


.80 


Fe xm 


203. 


.826 


27. 


.46 


25 


.02 


1, 


.10 


Fe xm 


204. 


.937 


7, 


.88 


8. 


.19 


0. 


.96 


Fexiv 


211. 


.316 


30. 


.71 


39 


.47 


0. 


.78 


Fexiv 


270. 


.519 


7. 


.97 


6. 


.96 


1. 


.14 


Fexiv 


274. 


.203 


17. 


.19 


18. 


.31 


0. 


.94 


Fexv 


284. 


.160 


28. 


.72 


21, 


.20 


1. 


.35 


Fexvi 


262. 


.984 


0. 


.17 


0. 


.42 


0. 


.42 



